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[i]  A  variety  of  spacebome  experiments  have  observed  polar  mesospheric  clouds  (PMC) 
since  the  late  20th  century.  Many  of  these  experiments  are  on  satellites  in  Sun- 
synchronous  orbits  and  therefore  allow  observations  only  at  fixed  local  times  (LT). 

Temperature  oscillations  over  the  diurnal  cycle  are  an  important  source  of  PMC  variability. 

In  order  to  quantify  long-term  natural  or  anthropogenic  changes  in  PMCs,  it  is  therefore 
essential  to  understand  their  variation  over  the  diurnal  cycle.  To  this  end,  we  employ  a 
prototype  global  numerical  weather  prediction  system  that  assimilates  satellite  temperature 
and  water  vapor  observations  from  the  ground  to  ~90  km  altitude.  We  assemble  the 
resulting  6  hourly  high-altitude  meteorological  assimilation  fields  from  June  2007  in  both 
LT  and  latitude  and  use  them  to  drive  a  one-dimensional  PMC  formation  model  with 
cosmic  smoke  serving  as  nucleation  sites.  We  find  that  there  is  a  migrating  diurnal 
temperature  tide  at  69°N  with  a  variation  of  ±4  K  at  83  km,  which  controls  the  variation 
of  PMC  total  ice  water  content  (IWC)  over  the  diurnal  cycle.  The  calculated  IWC  is 
normalized  to  observations  at  2300  LT  by  the  Solar  Occultation  for  Ice  Experiment  and 
allowed  to  vary  with  temperature  over  the  diurnal  cycle.  We  find  that  the  IWC  at  69°N 
has  a  single  maximum  between  0700  and  0800  LT  and  a  minimum  between  1900  and 
2200  LT  and  varies  by  at  least  a  factor  of  5.  The  calculated  variation  of  IWC  with  LT 
is  substantially  larger  at  57°N,  with  a  single  prominent  peak  near  0500  LT. 

Citation:  Stevens,  M.  H.,  et  al.  (2010),  Tidally  induced  variations  of  polar  mesospheric  cloud  altitudes  and  ice  water  content 
using  a  data  assimilation  system,  J.  Geophys.  Res .,  115 ,  D 18209,  doi:10.1029/2009JD0 13225. 


1.  Introduction 

[2]  Polar  mesospheric  clouds  (PMCs)  are  tenuous  layers 
of  ice  particles  that  form  near  83  km  altitude  over  the 
summer  polar  region.  When  observed  from  the  ground,  they 
are  also  known  as  noctilucent  clouds  (NLCs)  by  their  glowing 
appearance  as  observed  against  the  twilight  sky.  We  herein¬ 
after  refer  to  all  mesospheric  clouds  as  PMCs  for  consistency. 
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Some  analyses  of  satellite  observations  over  the  past  30  years 
suggest  that  PMCs  have  become  more  frequent  [e.g.,  Shettle 
et  al ,  2009]  and  brighter  [e.g.,  Thomas  et  al. ,  2003;  DeLand 
et  al ,  2003],  fueling  the  suggestion  that  they  are  indicators  of 
global  climate  change  [ Thomas  et  al,  1989;  Thomas ,  1996]. 

[3]  Ground-based  observations  indicate  that  PMCs 
exhibit  tidally  induced  variations  in  brightness,  altitude  and 
occurrence  frequency  [von  Zahn  etal.,  1998;  Chu  etal.,  2003, 
2006;  Fiedler  et  al,  2005].  On  the  other  hand,  PMC 
observations  are  typically  from  satellites  in  Sun- synchronous 
orbits  measuring  at  discrete  local  times  (LT)  during  any  given 
season.  This  is  important  because  ground-based  observations 
of  brightness  variations  at  a  fixed  location  vary  over  the 
diurnal  cycle  by  ±20%  or  more  [von  Zahn  et  al ,  1998;  Fiedler 
et  al,  2005],  yet  reported  multidecadal  PMC  albedo  trends 
are  less  than  1%/yr  [e.g.,  DeLand  et  al,  2007].  Constraining 
PMC  secular  variations  clearly  requires  a  quantitative  under¬ 
standing  of  diurnal  effects  so  that  inferred  trends  are  uncon¬ 
taminated  by  diurnal  variations  [Stevens  et  al ,  2007;  Kirkwood 
et  al,  2008].  This  is  especially  important  for  the  suite  of  solar 
backscattered  ultraviolet  (SBUV)  instruments,  which  have 
measured  PMCs  continuously  from  Sun- synchronous  orbits 
at  a  variety  of  fixed  LT  for  over  30  years.  For  Northern 
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Hemisphere  data,  these  SBUV  LT  are  biased  almost  exclu¬ 
sively  to  the  morning  from  about  1979  to  1990  and  then 
distributed  approximately  equally  over  the  diurnal  cycle  from 
1990  to  the  present  [ DeLand  et  al .,  2007]. 

[4]  One  useful  measure  of  PMC  variability  is  the  ice  mass 
column  density,  sometimes  called  the  total  ice  water  content 
(IWC)  [ Thomas  and  McKay ,  1985;  Stevens  et  al .,  2005; 
Rapp  and  Thomas ,  2006;  Stevens  et  al ,  2007;  Hervig  et  al , 
2009a,  2009b].  The  IWC  is  the  vertically  integrated  ice  mass 
density  and  therefore  a  direct  measure  of  the  total  ice  present 
in  a  PMC  layer.  A  quantitative  simulation  of  IWC  over  the 
diurnal  cycle  requires  knowledge  of  the  ambient  conditions 
with  LT. 

[5]  We  herein  determine  tidal  variations  of  temperature, 
water  vapor,  and  winds  over  the  diurnal  cycle  using  mete¬ 
orological  analyses  produced  by  the  data  assimilation 
component  of  the  Advanced-Level  Physics  High-Altitude 
(ALPHA)  Navy  Operational  Global  Atmospheric  Prediction 
System  (NOGAPS)  in  order  to  simulate  variations  in  IWC 
with  LT.  In  section  2,  we  describe  how  we  use  the  NOGAPS- 
ALPHA  analysis  for  June  2007  to  drive  a  one-dimensional 
microphysical  model  and  to  simulate  the  diurnal  variation  of 
IWC.  In  section  3,  we  compare  the  results  at  69°N  with  IWC 
observations  at  discrete  LT  by  the  Solar  Occultation  for  Ice 
Experiment  (SOFIE)  and  the  Cloud  Imaging  and  Particle  Size 
(CIPS)  instrument  on  the  NASA  Aeronomy  of  Ice  in  the 
Mesosphere  (AIM)  satellite  [Russell  et  al ,  2009;  Hervig  et 
al ,  2009a]  as  well  as  the  Student  Nitric  Oxide  Explorer 
(SNOE)  [Bailey  et  al ,  2005].  We  also  compare  the  model 
IWC  results  at  57°N  with  cloud  occurrence  frequency 
observations  over  the  diurnal  cycle  by  the  Spatial  Heterodyne 
Imager  for  Mesospheric  Radicals  (SHIMMER)  [Englert  et 
al ,  2008;  Stevens  et  al ,  2009].  In  section  4,  we  show  tidal 
variations  at  PMC  altitudes  over  three  successive  northern 
seasons  and  discuss  the  implications  of  all  our  results  for 
long-term  trend  analyses. 


2.  Modeling  Approach 

[6]  Our  modeling  approach  is  designed  to  quantify  the 
effects  of  migrating  tides  on  heavily  averaged  PMC  ob¬ 
servations  in  the  Northern  Hemisphere.  This  study  therefore 
does  not  address  the  geographical  variability  of  ice  layers 
[Berger  and  Liibken ,  2006],  hemispheric  differences  of 
PMCs  [e.g.,  Bailey  et  al ,  2007;  Liibken  and  Berger ,  2007]  or 
the  modeling  of  long-term  PMC  variability  [Liibken  et  al , 
2009].  We  focus  instead  on  how  PMCs  vary  with  LT  in  the 
Arctic  so  that  small  long-term  changes  inferred  from  existing 
satellite  data  sets  fixed  in  LT  may  be  uncoupled  from  diurnal 
variations.  This  may  benefit  the  identification  of  any  source 
of  long-term  PMC  change. 

[7]  We  first  describe  the  ambient  conditions  over  both 
latitude  and  LT  and  then  use  this  time  history  to  drive  a  one¬ 
dimensional  microphysical  model.  To  estimate  the  meridi¬ 
onal  transport  of  the  ice  particles,  we  use  back  trajectories 
based  upon  NOGAPS -ALPHA  analyzed  horizontal  winds. 
These  trajectories  and  hence  source  regions  vary  with  LT 
since  the  wind  field  varies  with  LT.  Starting  from  the  source 
region  for  the  target  LT  and  latitude  of  interest,  we  integrate 
forward  in  time  to  drive  the  one-dimensional  Community 
Aerosol  and  Radiation  Model  for  Atmospheres  (CARMA) 


[Rapp  and  Thomas ,  2006]  from  which  the  PMC  properties 
are  calculated. 

[8]  The  calculated  IWC  is  a  function  of  how  long  the  ice 
particles  are  subjected  to  temperatures  below  the  frost  point. 
This  particle  “lifetime”  may  be  related  to  processes  we  do 
not  consider,  such  as  gravity  wave  activity,  and  is  therefore 
not  well  constrained  by  our  model.  However,  we  use  the 
SOFIE  IWC  observations  as  a  point  of  reference  at  2300  LT, 
and  we  adjust  the  model  ice  particle  lifetime  to  agree  with 
these  data.  To  do  this  we  introduce  a  warm  bias  to  the  model 
and  cool  the  model  atmosphere  to  ambient  conditions  so  that 
the  resultant  ice  particle  growth  yields  an  IWC  equal  to 
SOFIE  IWC  observations  at  2300  LT.  To  model  the  diurnal 
variation  of  PMCs,  we  use  24  trajectories,  one  for  each  hour 
of  LT,  and  use  the  same  lifetime  for  each  simulation.  In  the 
following  two  subsections,  we  first  describe  our  composite 
wind  and  temperature  fields  from  NOGAPS-ALPHA.  We 
then  discuss  their  use  to  drive  CARMA. 

2.1.  NOGAPS-ALPHA 

[9]  NOGAPS-ALPHA  is  a  nonoperational  research  pro¬ 
totype  of  the  U.  S.  Navy’s  operational  global  numerical 
weather  prediction  system  that  extends  its  vertical  range  to 
~90  km  altitude,  allowing  us  to  simulate  the  PMCs  near 
83  km.  Unlike  climate  models,  this  system  is  focused  on 
accurate  “nowcasting”  and  short-term  forecasting  of  the 
global  atmosphere.  The  system  consists  of  two  main  com¬ 
ponents:  a  global  forecast  model  that  predicts  future  atmo¬ 
spheric  conditions  and  a  three-dimensional  variational  data 
assimilation  system  (NAVDAS:  NRL  Atmospheric  Varia¬ 
tional  Data  Assimilation  System)  that  provides  synoptic 
initial  conditions  for  those  forecasts  based  on  available  ob¬ 
servations.  The  two  components  work  synergistically  in  a 
coupled  forecast-assimilation  update  cycle,  with  forecasts 
providing  an  estimate  of  the  global  atmospheric  state  that  are 
updated  every  6  h  by  the  assimilation  of  global  observations. 
Our  focus  in  this  work  is  on  the  global  analysis  fields  in 
the  polar  summer  mesosphere  provided  by  the  assimilation 
of  mesospheric  observations  into  the  NOGAPS-ALPHA 
system. 

[10]  NOGAPS-ALPHA  analysis  fields  have  been  extended 
upward  by  assimilating  high-altitude  observations  from 
research  satellites,  in  addition  to  the  usual  complement  of 
operational  sensor  data  at  lower  altitudes.  Middle  atmo¬ 
sphere  temperature  and  water  vapor  observations  are  assim¬ 
ilated  from  the  Microwave  Limb  Sounder  (MLS)  on  the 
Aura  satellite  and  the  Sounding  of  the  Atmosphere  using 
Broadband  Emission  Radiometry  (SABER)  instrument  on 
the  Thermosphere  Ionosphere  Mesosphere  Energetics  and 
Dynamics  (TIMED)  satellite.  The  vertical  resolution  of  the 
SABER  temperature  observations  is  ~2  km  [Remsberg  et  al , 
2008]  and  the  vertical  resolution  of  the  MLS  temperature  and 
water  vapor  retrievals  is  ~  10-1 2  km  near  80  km  altitude 
[Lambert  et  al ,  2007;  Schwartz  et  al ,  2008].  Hoppel  et  al 
[2008]  described  initial  experiments  for  January  2006  that 
extended  the  analyses  to  0.01  hPa  by  assimilating  MLS  and 
SABER  temperatures.  Eckermann  et  al  [2009]  extended  the 
system  by  assimilating  newer  versions  of  the  temperature 
retrievals  up  to  0.002  hPa,  in  addition  to  MLS  water  vapor 
and  ozone.  It  is  these  latter  analysis  fields  that  we  use  in  the 
present  work. 
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[11]  The  effective  horizontal  resolution  of  this  NOGAPS- 
ALPHA  analysis  is  approximately  2°  latitude  x  2°  longitude, 
although  the  output  is  interpolated  onto  a  1°  x  1°  isobaric 
grid.  The  assimilated  temperature  fields  in  the  summer  polar 
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mesosphere  revealed  good  agreement  with  independent  sat¬ 
ellite  observations  [Eckermann  et  al .,  2009,  Figure  2].  Fur¬ 
thermore,  saturation  ratios  derived  diagnostically  from 
assimilated  temperature  and  humidity  fields  exhibited 
excellent  agreement  with  PMCs  observed  from  the  ground, 
from  SOFIE,  and  from  SHIMMER  [Eckermann  etal.,  2009]. 

[12]  The  geophysical  variability  of  the  horizontal  winds  is 
~20  m/s  from  one  6  h  update  cycle  to  the  next  and  the 
geophysical  variability  of  the  temperatures  is  ~5  K,  which 
combine  to  yield  unreasonably  large  uncertainties  in  our 
back  trajectories  and  cloud  simulations.  To  reduce  the  var¬ 
iability,  we  average  the  NOGAPS-ALPHA  synoptic  fields 
of  temperature,  water  vapor,  and  horizontal  and  vertical 
winds  on  a  LT  grid  for  all  of  June  2007.  This  month  is 
chosen  because  it  is  during  the  TIMED  yaw  cycle  that  in¬ 
cludes  north-looking  SABER  data,  thereby  maximizing  the 
amount  of  assimilated  temperature  data  in  the  polar  summer. 
This  June  average  of  120  6  hourly  assimilation  update 
cycles  removes  all  longitudinal  variability  and  instead  yields 
fields  that  isolate  the  migrating  tides,  which  are  now 
described  as  a  function  of  latitude  and  LT  only.  In  this  way, 
we  derive  representative  conditions  for  the  polar  summer 
mesosphere  for  use  in  comparisons  with  longitudinally  and 
seasonally  averaged  PMC  observations. 

[13]  The  average  diurnal  variations  for  June  2007  from 
NOGAPS-ALPHA  are  shown  at  83  km  geometric  altitude 
in  Figure  la  for  temperature  and  Figure  lb  for  water  vapor 
with  the  calculated  degree  of  supersaturation  shown  in 
Figure  lc.  Although  the  equilibrium  vapor  pressure  of 
water  vapor  over  ice  at  such  low  temperatures  has  not  been 
measured  directly,  our  diagnostic  supersaturation  ratios  in 
Figure  lc  use  the  expression  of  Marti  and  Mauersberger 
[1993]  as  recommended  by  Rapp  and  Thomas  [2006]. 
From  Figure  la,  the  temperature  at  69°N  is  on  average 
143  K,  in  good  agreement  with  the  climatology  of  Lubken 
[1999],  who  finds  the  temperature  to  be  146  ±  4  K  at  this 
altitude  in  mid- June.  A  predominantly  diurnal  oscillation 
is  present  in  temperature  with  an  amplitude  of  about  4  K 
near  69 °N,  a  minimum  near  0300  LT  and  a  maximum  near 
1700  LT,  in  agreement  with  the  analysis  at  60°N  of 
Eckermann  et  al  [2009,  Figure  13].  This  result  from 
NOGAPS-ALPHA  is  also  in  good  agreement  with  inde¬ 
pendent  ground-based  temperature  observations  by  Singer 
et  al  [2003],  who  report  that  the  diurnal  tide  dominates 
summer  mesospheric  temperature  variability  at  Arctic  lati¬ 
tudes  with  an  amplitude  of  3-8  K.  Near  the  equator  at  82  km 
altitude,  temperature  tides  from  NOGAPS-ALPHA  have 
larger  amplitudes  of  5-10  K,  in  general  agreement  with 
data-validated  results  from  the  Canadian  Middle  Atmo¬ 
sphere  Model  for  June  [ McClandress ,  2002;  Eckermann  et 
al,  2009]. 


Figure  1.  NOGAPS-ALPHA  analyses  at  83  km  geometric 
altitude  and  for  June  2007.  The  results  for  (a)  temperature, 
(b)  water  vapor,  and  (c)  degree  of  supersaturation  (N)  are 
averaged  in  LT  from  50°N-85°N.  Latitudes  poleward  of 
85°N  are  omitted  due  to  the  lack  of  satellite  data.  Local  mid¬ 
night  is  at  the  bottom  of  each  plot  and  noon  is  at  the  top  so 
that  LT  progresses  counterclockwise  around  each  image. 
Note  that  PMCs  can  exist  where  S  >  1. 
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[14]  The  water  vapor  in  Figure  lb  shows  a  slight  overall 
enhancement  near  65 °N  relative  to  mean  values  at  other 
latitudes.  This  enhancement  near  65 °N  was  predicted  in  a 
modeling  study  by  von  Zahn  and  Berger  [2003]  as  being  due 
to  equatorward  transport  of  PMCs  to  warmer  regions  where 
the  ice  sublimates;  however,  they  were  unable  to  verify  this 
prediction  with  observations.  Since  the  NOGAPS-ALPHA 
synoptic  water  vapor  analysis  is  driven  by  MLS  data,  our 
result  in  Figure  lb  can  be  considered  as  a  validation  of  the  von 
Zahn  and  Berger  prediction. 

[15]  The  water  vapor  mixing  ratio  from  NOGAPS- 
ALPHA  in  Figure  lb  is  about  5  ppmv  at  69°N  and  2300  LT. 
This  is  in  good  agreement  with  independent  SOFIE  ob¬ 
servations  fixed  at  2300  LT  that  show  that  the  water  vapor 
mixing  ratio  is  about  5  ppmv  in  June  2007  near  83  km  altitude 
[ Gordley  et  al. ,  2009].  The  variation  of  the  NOGAPS- 
ALPHA  water  vapor  over  the  diurnal  cycle  is  about 
±0.25  ppmv  in  Figure  lb.  This  amplitude  is  smaller  than 
that  modeled  by  von  Zahn  and  Berger  who  show  variations 
of  about  ±3  ppmv  at  82  km.  This  difference  may  be  due  to 
the  coarse  vertical  resolution  (10-12  km)  of  the  MLS  water 
vapor  profiles  used  in  the  data  assimilation. 

[16]  The  degree  of  supersaturation  (S)  calculated  from  the 
temperature  and  water  vapor  in  Figures  la  and  lb,  respec¬ 
tively,  is  shown  in  Figure  lc.  PMCs  can  exist  in  regions 
where  S  >  1  and  the  low  morning  temperatures  in  Figure  la 
are  reflected  in  the  high  degree  of  saturation  during  the  same 
time  period  in  Figure  lc.  Note  that  S  >  1  at  all  LT  poleward 
of  ~60°N. 

[  1 7]  The  derived  zonal,  meridional,  and  vertical  winds  from 
the  NOGAPS-ALPHA  analysis  are  shown  in  Figures  2a-2c. 
The  vertical  winds  are  calculated  diagnostically  from  the 
divergence  of  the  horizontal  wind  fields  [Hogan  and 
Rosmond ,  1991]  and  peak  near  0400  and  1500  LT.  The 
importance  of  vertical  winds  in  the  calculation  of  meso¬ 
spheric  ice  particle  properties  is  discussed  in  detail  by  Berger 
and  von  Zahn  [2002].  Note  that  the  NOGAPS-ALPHA 
absolute  vertical  wind  speeds  reach  up  to  20  cm/s.  This  is 
considerably  larger  than  those  used  by  some  other  previous 
PMC  modeling  studies,  which  show  mean  vertical  winds 
from  2  to  5  cm/s  near  83  km  altitude  [e.g.,  von  Zahn  and 
Berger ,  2003;  Siskind  et  al .,  2005;  Herbort  et  al.,  2007]. 
We  note  that  ground-based  lidar  observations  near  40°N 
have  indicated  a  predominantly  semidiurnal  variation  in 
vertical  winds  with  amplitudes  near  30  cm/s  in  this  region  of 
the  atmosphere  [Kwon  et  al. ,  1987].  For  comparison, 
average  sedimentation  rates  of  PMC  ice  particles  at  83  km 
altitude  are  smaller  and  around  0.1-6  cm/s  [. Klostermeyer , 
1998]. 

[is]  To  distinguish  the  relative  amplitudes  of  the  migrating 
diurnal  and  semidiurnal  tidal  component,  we  have  Fourier 
transformed  the  temperatures  shown  in  Figure  la  and  show 
the  results  in  Figures  3a  and  3b.  For  these  results  we  have 


Figure  2.  NOGAPS-ALPHA  analysis  results  for  (a)  zonal 
winds,  (b)  meridional  winds,  and  (c)  derived  vertical  winds 
averaged  in  LT  over  June  2007  at  83  km  geometric  altitude. 
Positive  zonal  winds  are  eastward,  positive  meridional  winds 
are  northward,  and  positive  vertical  winds  are  upward.  The 
zero  wind  contours  are  overplotted  in  Figures  2b  and  2c. 
Latitudes  poleward  of  85°N  are  omitted  due  to  lack  of  data. 
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Figure  3.  Mean  variation  with  local  time  and  latitude  of  temperature  at  83  km  geometric  altitude  due  to 
(a)  the  migrating  diurnal  tide  and  (b)  the  migrating  semidiurnal  tide.  Results  are  monthly  means  for  June 
2007.  See  text  for  details. 


assumed  that  a  24  h  period  corresponds  exclusively  to  the 
wave  number  1  solar  migrating  diurnal  tide  (Figure  3  a)  and 
that  a  12  h  period  corresponds  exclusively  to  the  wave 
number  2  solar  migrating  semidiurnal  tide  (Figure  3b). 
Figure  3  a  shows  that  the  diurnal  tide  has  the  largest 
amplitude  between  50°N  and  60 °N  and  dominates  over  the 
semidiurnal  tide  in  Figure  3b  throughout  the  PMC  region. 

[19]  To  illustrate  how  pervasive  the  oscillation  of 
NOGAPS-ALPHA  vertical  winds  is  with  altitude  over  the 
diurnal  cycle,  Figure  4a  shows  a  cross  section  of  the  vertical 
winds  at  69°N  over  both  LT  and  altitude  in  the  PMC  altitude 
region  between  80  and  87  km.  The  calculated  PMC  altitudes 
using  CARMA  are  also  shown  and  are  determined  by  the 
altitude  of  maximum  ice  volume  density  for  each  simu¬ 
lation.  The  PMC  altitude  varies  over  the  diurnal  cycle  but 
on  average  is  near  83  km,  consistent  with  observations  and 
model  results  [Liibken  et  al .,  2008].  As  expected,  the  PMC 
altitudes  reach  a  maximum  at  the  node  between  the  max¬ 
imum  upward  and  downward  vertical  wind.  Figure  4a  also 
illustrates  that  the  strong  upward  winds  near  0400  and 
1500  LT  are  present  throughout  the  upper  mesosphere 
and  would  help  drive  the  altitude  variations  of  both  PMC 
and  polar  mesospheric  summer  echoes  (PMSE)  in  a  similar 
way  as  is  observed  from  ground-based  observations  [ von 
Zahn  and  Bremer ,  1999;  Hoffmann  et  al .,  2008].  Direct 
comparisons  of  our  model  results  in  Figure  4  to  altitude 
variations  over  the  diurnal  cycle  from  ground-based  observa¬ 
tions  are  limited  by  differences  in  the  ice  particle  history 
inherently  present  in  observations  at  a  fixed  location  such 
as  gravity  wave  activity,  planetary  wave  activity,  and  non¬ 
migrating  tides. 

[20]  Figure  4b  shows  the  variations  in  temperatures  over 
the  diurnal  cycle  for  the  same  altitude  region.  Note  that  the 
PMC  altitude  variations  driven  by  the  vertical  winds  slightly 
modify  the  diurnal  oscillation  of  temperatures  on  the  ice 
particles  compared  to  oscillations  at  a  fixed  altitude  surface 
of  83  km.  So  although  the  temperature  excursions  at  83  km 
are  ±4  K,  the  temperature  excursions  to  which  ice  particles 
are  subjected  are  about  ±3  K. 

[21]  In  order  to  help  validate  the  horizontal  and  vertical 
winds  shown  in  Figure  4a,  we  compare  them  with  ground 
based  horizontal  meteor  winds  observed  at  69°N  (Andenes, 
Norway)  and  85  km  geometric  altitude  in  Figure  5.  In 
addition,  we  performed  five  short  term  integrations  of  the 
forecast  model  component  of  NOGAPS-ALPHA,  initialized 


at  different  days  spanning  the  month  of  June  2007,  and 
present  the  calculated  winds  from  the  third  model  day  of 
those  simulations  as  an  additional  curve.  A  semidiurnal 
oscillation  in  both  the  meridional  and  zonal  component, 
with  a  stronger  peak  later  in  the  day,  is  evident  in  all  three 
curves.  The  agreement  between  both  components  of 
NOGAPS-ALPHA  (forecast  model  and  assimilations)  and 
with  the  observations  is  very  satisfying  because  the  two 
NOGAPS-ALPHA  components  derive  winds  differently. 
NOGAPS-ALPHA  does  not  assimilate  any  middle  atmo¬ 
spheric  wind  data  directly,  but  rather,  as  part  of  the  tem¬ 
perature  assimilation,  calculates  correlated  temperature  and 
wind  increments  based  upon  a  gradient  wind  approximation 
in  the  off-diagonal  elements  of  the  background  error  covari¬ 
ance  matrix.  The  forecast  model  winds  are  constrained  by  the 
physical  parameterizations  in  the  model  (e.g.,  gravity  wave 
drag,  diffusion  and  tides  forced  by  the  model  at  lower 
altitudes  that  propagate  into  the  mesosphere).  Thus  the 
NOGAPS-ALPHA  winds  are  primarily  data  driven;  the 
forecast  winds  are  primarily  constrained  by  the  model 
physics.  The  validation  presented  in  Figure  5  with  the 
meteor  radar  winds  can  thus  be  considered  a  true  three-way 
intercomparison  between  two  approaches  to  calculating 
winds  and  the  Andenes  data. 

2.2.  Ice  Particle  Trajectories  and  CARMA 

[22]  As  noted  earlier,  our  approach  is  a  two-step  process. 
We  first  use  the  NOGAPS-ALPHA  analyzed  horizontal 
winds  to  calculate  hypothetical  cloud  parcel  trajectories  at 
83  km  altitude  from  a  given  latitude.  We  focus  herein  on 
two  latitudes  for  which  we  have  data  in  June  2007:  69°N 
and  57°N.  Contributions  from  the  effects  of  vertical  winds 
are  small  so  variations  on  the  trajectory  due  to  changes  in 
the  ice  particle  altitude  are  not  included  here.  Since  the 
prevailing  motion  of  the  ice  particle  is  from  east  to  west,  a 
small  correction  is  made  to  the  LT  at  each  time  step  because 
the  zonal  winds  are  pushing  the  ice  parcels  toward  earlier 
LT.  For  an  average  westward  wind  of  30  m/s,  this  correction 
amounts  to  about  10  min/h  at  69 °N. 

[23]  The  second  step  uses  CARMA  to  simulate  PMC 
formation  and  growth.  We  integrate  forward  in  time  along 
the  trajectory  and  use  the  appropriate  temperature  and  ver¬ 
tical  wind  fields  at  each  time  step  as  inputs  to  CARMA. 
CARMA  has  been  used  in  many  previous  PMC  studies,  and 
we  use  the  same  version  as  that  used  in  the  work  of  Rapp 
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Figure  4.  (a)  Time-height  cross  section  of  vertical  winds 
(contour  labels  in  cm/s)  in  the  upper  mesosphere  from  the 
NOGAPS-ALPHA  analysis  averaged  over  June  2007  at 
69°N.  The  winds  are  directed  both  upward  and  downward 
with  the  strongest  upwelling  near  0300  and  1500  LT,  as  indi¬ 
cated  in  Figure  2c.  Calculated  PMC  altitudes  are  shown, 
indicating  the  strong  correlation  with  vertical  winds,  (b)  Same 
as  in  Figure  4a,  but  for  temperature  (contour  labels  in  K). 

and  Thomas  [2006],  although  modified  to  use  the  time- 
dependent  inputs  along  its  trajectory  as  described  above.  We 
use  the  same  eddy  diffusion  as  recommended  by  Rapp  and 
Thomas  and  do  not  vary  it  with  LT.  Ice  properties  are  saved 
from  the  last  time  step  at  each  hour  of  LT  from  a  fixed 
(target)  latitude  for  comparison  to  observations  over  the 
diurnal  cycle  at  that  latitude. 

[24]  The  NOGAPS-ALPHA  temperature,  water  vapor, 
and  vertical  winds  are  provided  in  1 0  longitude  increments, 
equivalent  to  every  4  min  in  LT.  Although  the  microphysics 
is  recalculated  every  100  s,  the  time  increment  over  which 
CARMA  is  updated  with  new  temperature  and  vertical  wind 
profiles  is  therefore  every  240  s.  The  water  vapor  is  initialized 
once,  using  the  relevant  profile  for  the  target  local  time  and 
latitude,  and  allowed  to  redistribute  so  that  the  water  vapor 


responds  to  vertical  transport  rather  than  horizontal  transport. 
The  vertical  grid  spacing  for  NOGAPS-ALPHA  is  about 
2  km  and  this  output  is  interpolated  onto  the  finer  (0.25  km) 
grid  of  CARMA. 

[25]  We  assume  that  cosmic  smoke  is  the  nucleation 
source  of  all  the  ice  particles  and  initialize  the  model  once 
with  the  distribution  reported  for  mid- June  at  74°N  by 
Bardeen  et  al.  [2010].  These  model  smoke  results  show  that 
it  is  depleted  in  the  summer  polar  mesosphere  due  to  rela¬ 
tively  rapid  equatorward  transport  [ Megner  et  al .,  2008; 
Bardeen  et  al .,  2008],  in  agreement  with  SOFIE  smoke 
observations  [Hervig  et  al .,  2009c]. 

[26]  We  have  chosen  to  investigate  PMC  IWC  since  it  is  a 
vertically  integrated  property,  thereby  avoiding  the  need  to 
accurately  simulate  the  peak  of  a  narrow  ice  layer,  and  since 
it  is  relatively  insensitive  to  variations  in  the  ambient  con¬ 
ditions  other  than  temperature  and  water  vapor  [Rapp  and 
Thomas ,  2006,  Table  1].  We  first  present  a  sensitivity  study 
showing  how  IWC  varies  depending  on  the  modeling 
approach.  Specifically,  we  show  how  IWC  varies  with  fixed 
ambient  conditions,  with  variations  over  latitude  and  LT  and 
with  varying  ice  particle  lifetimes.  These  simulations  focus 
on  results  near  2300  LT  since  this  corresponds  to  the  SOFIE 
observations.  The  SOFIE  IWC  is  an  average  for  June  2007 
(39  yug/m2),  but  we  increase  this  by  39%  (54  yug/m2)  based 
on  the  reported  latitudinal  variation  of  ice  mass  from  the 
average  SOFIE  latitude  (66°N)  in  June  to  69°N  [Stevens  etal ., 
2007]. 

[27]  Figure  6  shows  a  variety  of  CARMA  simulations  in 
which  we  have  varied  the  assumptions  used  in  the  calcula¬ 
tion  of  IWC.  Each  of  the  simulations  in  Figure  6  is  inte¬ 
grated  for  48  h,  but  only  results  for  the  last  24  h  are  shown. 
Because  the  IWC  is  a  function  of  the  ice  particle  growth  rate 
and  the  particle  lifetime,  we  introduce  a  warm  bias  to  the 
ambient  temperatures  and  remove  the  bias  in  different  ways 
as  described  below.  The  water  vapor  profile  for  all  simula¬ 
tions  in  Figure  6  is  initialized  once  and  allowed  to  redis¬ 
tribute  vertically  through  freeze  drying. 


Figure  5.  A  comparison  of  (top)  meridional  and  (bottom) 
zonal  winds  from  the  NOGAPS-ALPHA  analysis  and  the 
NOGAPS-ALPHA  forecast  with  ground  based  meteor  wind 
observations  during  the  same  June  2007  time  period  at  the 
geometric  altitude  of  85  km. 
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[28]  Figure  6a  shows  an  IWC  simulation  where  the 
ambient  temperature  is  given  a  30  K  warm  bias  and  then 
abruptly  set  to  the  temperature  profile  at  69°N  and  2300  LT, 
after  which  the  ice  particles  are  allowed  to  nucleate  and 
grow  over  a  period  of  24  h.  This  simple  approach  also  uses  a 


2.8 


71.4 


Locol  Time  (h) 

7.5  1 2.3  1 7.5 


Lotitude  (N) 
70.3  70.3 


constant  vertical  wind  profile  from  69 °N  and  therefore 
neglects  all  effects  of  diurnal  variations  and  horizontal  trans¬ 
port  from  higher  latitudes.  For  the  simulation  in  Figure  6a, 
we  calculate  an  IWC  of  140  yug/m2  at  the  end  of  the  final 
24  h.  Rapp  and  Thomas  [2006]  use  the  same  approach  of  a 
“cold  start”  with  a  fixed  temperature  and  vertical  wind 
profile  to  calculate  the  IWC.  Although  they  did  not  use  the 
same  temperatures,  vertical  winds,  water  vapor,  or  meteoric 
smoke  as  our  simulation,  they  nonetheless  find  an  IWC  of 
135  yug/m2  for  their  reference  case  after  24  h,  very  similar 
to  our  calculated  value.  These  results  are  however  a  factor 
of  ~2.5  larger  than  the  SOFIE  observation  indicated  in 
Figure  6a.  Bardeen  et  al.  [2010]  used  the  three-dimensional 
Whole-Atmosphere  Community  Climate  Model  (WACCM) 
with  microphysics  from  CARMA  to  derive  an  IWC  of  12- 
51  //g/m2  at  66°N  and  2300  LT,  in  good  agreement  with  the 
SOFIE  observations  for  June  2007. 

[29]  Figure  6b  shows  the  results  of  two  simulations,  each 
of  which  includes  the  effects  of  transport  from  higher  lati¬ 
tudes.  The  variations  of  LT  and  latitude  over  one  day  of 
elapsed  time  are  indicated  on  the  upper  axes  of  this  panel.  In 
one  case,  we  have  started  the  ice  formation  abruptly  with  the 
temperature  and  vertical  wind  profiles  starting  24  h  prior 
and  run  the  simulation  for  24  h  while  continuously  updating 
these  profiles  along  the  trajectory.  This  yields  an  IWC  that  is 
about  a  factor  of  4  larger  than  the  SOFIE  observation. 

[30]  This  “cold  start”  approach,  whereby  CARMA  is 
initialized  from  a  highly  supersaturated  state  and  with 
clouds  growing  quickly,  is  likely  unrealistic.  It  is  more 
likely  that  clouds  grow  more  slowly  as  the  air  starts  out 
subsaturated,  then  cools,  and  becomes  progressively  more 
saturated.  Unfortunately  we  cannot  fully  address  this  since 
our  microphysical  model  is  not  coupled  to  the  background 
atmosphere  (as  in  the  work  of  Bardeen  et  al.  [2010]).  Fur¬ 
thermore,  there  is  little  observational  data  as  to  the  time 
history  of  the  PMCs.  Instead,  our  approach  is  to  simply 
assume  the  air  is  subsaturated  and  apply  a  progressively 
decreasing  temperature  bias  along  the  parcel  trajectory.  To 
fit  the  SOFIE  observations,  we  set  this  cooling  rate  to 


Figure  6.  (a)  The  calculated  time  variation  of  IWC. 
NOGAPS -ALPHA  temperature,  water  vapor,  and  vertical 
wind  profiles  are  used  for  69°N  and  2300  LT.  The  temper¬ 
ature  and  vertical  wind  profiles  are  fixed,  (b)  Two  simula¬ 
tions  showing  the  variation  of  IWC  with  latitude  and  LT 
for  an  ice  particle  advected  through  ambient  temperatures 
and  vertical  winds  to  69°N.  The  dashed  line  shows  a  24  h 
simulation  where  the  temperature  profiles  from  NOGAPS- 
ALPHA  are  used  throughout  (“Cold  Start”).  The  solid  line 
introduces  a  temperature  bias  that  is  gradually  removed 
(see  text),  thereby  reducing  the  ice  particle  growth  and  the 
IWC  (“Short  Lifetime”).  The  dotted  portion  of  the  Short 
Lifetime  case  indicates  the  time  period  during  which  the 
ambient  conditions  have  not  yet  been  reached,  (c)  Similar 
to  the  Short  Lifetime  case  of  Figure  6b  except  three  other 
LTs  are  added  for  comparison.  The  indicated  LTs  corre¬ 
spond  to  the  last  time  step  of  each  simulation,  (d)  The 
IWC  variation  at  2300  and  0700  LT  using  two  different  cool¬ 
ing  rates  from  an  initial  temperature  bias:  -0.5  K/h  (22.8  K 
bias)  and  -1.0  K/h  (44.4  K  bias).  The  temperature  biases 
are  selected  so  that  the  SOFIE  data  at  2300  LT  is  fit. 
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Figure  7.  The  calculated  IWC  over  the  diurnal  cycle  at 
69°N  for  June  2007.  Three  different  solutions  are  pre¬ 
sented  representing  three  assumed  parcel  cooling  rates:  slow 
(0.50  K/h),  medium  (0.68  K/h),  and  fast  (1.00  K/h),  where 
the  shading  represents  the  range  within  these  solutions 
(see  text).  AIM/CIPS  (blue)  and  AIM/SOFIE  (red)  data 
from  the  same  time  period  are  overplotted.  The  CIPS  obser¬ 
vations  are  averaged  over  a  3°  latitude  band  around  69°N 
and  have  been  scaled  up  to  be  consistent  with  simultaneous 
SOFIE  observations  in  separate  common  volume  observa¬ 
tions.  Also  shown  are  observations  from  SNOE,  which  were 
collected  between  0900  and  1100  LT  and  have  been  scaled 
upwards  for  representative  solar  minimum  conditions. 

-0.68  K/h  from  an  initial  temperature  bias  of  +30  K  at  all 
altitudes  in  Figure  6b.  The  ice  particles  are  therefore  sub¬ 
jected  to  ambient  conditions  for  only  3—4  h  out  of  a  total  of 
48  h  of  microphysics.  We  note  that  this  effective  ice  particle 
lifetime  of  3-4  h  is  an  order  of  magnitude  smaller  than  the 
ice  particle  lifetime  inferred  by  Berger  and  von  Zahn  [2007], 
who  instead  derived  their  lifetime  from  the  observed  size 
distribution  at  the  peak  of  the  ice  layer  at  69°N.  However, 
the  particle  lifetime  used  herein  is  consistent  with  the  more 
recent  study  of  Zasetsky  et  al.  [2009],  who  argued  that 
particles  form  in  2-20  h.  As  indicated  in  Figure  6b,  longer 
lifetimes  increase  the  IWC  beyond  the  SOFIE  observations 
for  our  June  2007  case  study  at  2300  LT. 

[31]  Figure  6c  shows  the  results  for  simulations  at  four 
different  LTs:  0700  LT,  1400  LT,  2100  LT,  and  2300  LT. 
0700  LT  and  2100  LT  are  chosen  because  they  yield  the 
maximum  and  minimum  IWC  over  the  diurnal  cycle,  respec¬ 
tively,  1400  LT  is  chosen  because  it  corresponds  to  some  of 
the  CIPS  observations  (discussed  in  section  3)  and  2300  LT 
corresponds  to  the  SOFIE  observations  reproduced  from 
Figure  6b.  Figure  6c  illustrates  that  the  simulations  that  end 
at  0700  and  1400  LT  produce  clouds  earlier  than  the  others, 
but  that  the  ambient  conditions  are  reached  at  the  same  point 
in  each  simulation.  The  IWC  is  largest  at  0700  LT  because 
this  is  just  after  the  low  temperature  phase  of  the  diurnal 
cycle  and  IWC  is  smallest  at  2100  LT  because  this  is  just 
after  the  time  of  highest  temperatures  (see  Figure  4). 

[32]  The  cooling  rate  (-0.68  K/h)  and  the  temperature  bias 
(+30  K)  used  in  Figure  6c  are  selected  to  yield  a  solution 
consistent  with  SOFIE  at  2300  LT.  However,  other  combi¬ 


nations  of  these  two  parameters  can  also  fit  the  SOFIE  data. 
We  herein  also  consider  cooling  rates  of -0.5  and  -1.0  K/h  to 
assess  the  impact  on  the  resultant  IWC,  consistent  with 
representative  rates  for  heterogeneous  nucleation  as  dis¬ 
cussed  in  the  recent  work  of  Murray  and  Jensen  [2010]  that 
bracket  our  baseline  case.  These  cooling  rate  changes  vary 
the  exposure  time  to  conditions  which  favor  growth.  The 
temperature  biases  for  these  two  cases  are  chosen  to  fit  the 
SOFIE  data  at  2300  LT  at  the  end  of  each  48  h  integration. 
Figure  6d  shows  the  results  for  a  -0.5  K/h  cooling  rate 
(+22.8  K  bias)  and  -1.0  K/h  cooling  rate  (+44.4  K  bias)  at 
2300  LT  and  also  at  the  maximum  IWC  at  0700  LT.  The 
IWC  at  0700  LT  are  very  similar,  and  we  consider  the 
effects  of  these  different  cooling  rates  over  the  entire  diurnal 
cycle  in  section  3. 

3.  Diurnal  IWC  Simulations 
3.1.  Results  for  69°N 

[33]  Figure  7  shows  the  calculated  IWC  at  each  hour  of 
LT  over  the  diurnal  cycle.  For  each  of  the  24  simulations, 
the  temperatures  and  the  vertical  winds  were  specified  along 
a  trajectory  defined  by  the  horizontal  winds  and  ending  at 
the  local  time  of  interest.  The  water  vapor  is  initialized  once 
for  each  LT  at  69 °N  and  allowed  to  redistribute  vertically 
through  freeze  drying,  as  for  the  simulations  in  Figure  6.  In 
addition  to  our  baseline  case  shown  in  Figure  6c,  we  include 
results  with  more  cooling  and  less  cooling  as  discussed  in 
section  2.  The  shaded  area  in  Figure  7  shows  the  range  of 
solutions  at  each  LT  from  the  three  different  simulations.  As 
with  the  PMC  altitudes  shown  in  Figure  4,  a  comparison  of 
the  IWC  variation  over  the  diurnal  cycle  with  ground-based 
observations  from  a  single  location  is  misleading  because 
the  ice  particles  observed  at  one  location  can  be  subject  to 
processes  that  we  cannot  reproduce  with  our  approach.  We 
instead  compare  the  simulated  IWC  to  zonally  averaged 
satellite  observations  measuring  at  discrete  LT. 

[34]  On  Figure  7  there  are  also  two  zonally  averaged 
observations  of  IWC  from  CIPS  (version  3.20  retrievals) 
[see  Rusch  et  al ,  2009;  Benze  et  al ,  2009],  which  measures 
near  1350  and  2240  LT  in  the  Northern  Hemisphere  on  the 
descending  and  ascending  nodes  of  the  AIM  orbit,  respec¬ 
tively.  SOFIE  is  a  solar  occultation  instrument  and  sensi¬ 
tively  measures  the  IWC  directly  whereas  CIPS  is  a  nadir 
imager  and  less  sensitive,  so  we  have  applied  an  upward 
adjustment  to  the  CIPS  data  using  separate  common  volume 
(CV)  observations  [Russell  et  al ,  2009].  This  upward 
adjustment  is  discussed  further  below. 

[35]  Near  69°N,  SOFIE  and  CIPS  execute  observations 
that  are  specifically  coordinated  to  observe  the  same  volume 
of  ice.  We  define  the  average  IWC  from  these  CV  observa¬ 
tions  as  the  vertical  integral  of  the  ice  mass  density  through 
a  cloud,  multiplied  by  the  cloud  frequency  during  June 
2007.  In  the  case  of  CIPS,  the  cloud  frequency  is  determined 
by  the  number  of  CIPS  pixels  in  which  a  cloud  was  detected 
divided  by  the  total  number  of  CIPS  pixels  along  the  SOFIE 
line  of  sight.  For  SOFIE,  the  cloud  frequency  is  the  number 
of  occultations  in  which  a  cloud  was  detected  divided  by  the 
total  number  of  occultations.  Using  this  approach,  SOFIE 
yields  an  IWC  consistently  higher  than  CIPS  because  many 
clouds  along  the  SOFIE  line  of  sight  are  below  the  CIPS 
detection  threshold. 
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[36]  To  account  for  the  reduced  sensitivity  of  CIPS,  we 
scale  the  CIPS  CV  IWC  up  to  the  SOFIE  CV  IWC.  How¬ 
ever,  we  do  not  know  independently  how  much  IWC  is 
below  the  CIPS  threshold  away  from  the  CV.  We  therefore 
calculate  a  scale  factor  for  the  undetectable  ice  in  two  ways: 
one  scale  factor  (4.9)  assumes  that  the  observations  with  no 
PMCs  are  indeed  clear  air  and  one  (1.9)  assumes  that  these 
observations  are  at  the  CIPS  threshold.  We  find  that  the 
average  CIPS  IWC  threshold  is  18  /xg/m2,  which  corre¬ 
sponds  on  average  to  a  CIPS  directional  albedo  threshold  of 
5  x  10~6  sr_1.  To  be  consistent  with  this  approach  in  the  CV, 
at  69°N  we  assign  values  of  zero  where  there  is  clear  air  and 
use  the  scale  factor  of  4.9  for  one  solution  and  the  smaller 
scale  factor  of  1 .9  for  the  other  solution.  The  results  are  shown 
in  Figure  7,  and  on  average  the  CIPS  IWC  near  1400  LT  is 
58  /ig/m2  whereas  the  IWC  near  2300  LT  is  31%  smaller  at 
40  yug/m2.  The  agreement  of  the  scaled  CIPS  data  with  the 
simulations  in  Figure  7  is  good,  although  these  data  do  not 
significantly  help  to  validate  the  large  dynamic  range  of 
the  IWC  simulated  by  CARMA  using  NOGAPS-ALPHA 
analyses  because  CIPS  does  not  make  measurements  at 
those  LTs  when  our  model  predicts  a  large  IWC. 

[37]  To  help  validate  the  larger  morning  IWC,  we  use 
observations  from  SNOE  near  1000  LT.  SNOE  was  launched 
into  a  Sun- synchronous  orbit  in  February  1998  and  observed 
PMCs  for  nearly  six  years  over  both  poles  [ Bailey  et  al, 
2005].  Stevens  et  al  [2007]  inferred  the  PMC  ice  mass 
from  SNOE  observations  in  the  Arctic  from  the  cloud 
brightnesses  and  occurrence  frequencies.  The  reported 
SNOE  ice  mass  for  midsolar  cycle  conditions  at  70°N  is 
equivalent  to  an  IWC  of  53-106  yug/m2,  depending  on  the 
ice  particle  size  distribution.  Since  the  study  here  focuses  on 
solar  minimum  conditions  and  since  PMCs  are  known  to  be 
brighter  and  more  frequent  at  solar  minimum,  we  apply  an 
upward  correction  to  this  IWC,  which  is  described  below. 

[38]  Several  studies  have  calculated  the  solar  cycle  vari¬ 
ation  of  ice  content  (sometimes  called  H20(ice)),  ice  mass  or 
IWC,  which  all  scale  the  same  way.  Siskind  et  al  [2005] 
reported  a  ratio  of  the  ice  content  from  solar  minimum  to 
maximum  of  a  factor  of  3  near  70°N,  which  they  regarded  as 
a  likely  upper  limit.  Stevens  et  al  [2007]  used  SBUV  data  to 
show  that  the  variation  of  ice  mass  is  at  least  a  factor  of 
2  over  the  solar  cycle.  Using  these  studies,  we  therefore  take 
the  solar  cycle  variation  of  IWC  to  range  from  a  factor  of 
2-3  from  solar  maximum  to  solar  minimum.  We  note  that 
this  is  somewhat  greater  than  the  IWC  variation  over  the 
solar  cycle  inferred  in  the  modeling  study  of  Bardeen  et  al 
[2010],  who  determined  that  at  solar  minimum  the  IWC  is 
1 .44  times  greater  than  at  solar  maximum.  To  conservatively 
represent  the  range  of  possible  solutions,  we  therefore  apply 
a  factor  of  1 .4  for  our  lower  IWC  limit  from  midsolar  cycle 
to  solar  minimum  and  a  factor  of  2.1  for  our  upper  IWC 
limit  to  solar  minimum.  We  find  that  the  scaled  SNOE 
observations  range  from  74  to  223  yug/m2.  This  is  the  range 
overplotted  in  Figure  7  and  is  in  agreement  with  the  IWC 
calculation.  The  calculated  difference  between  the  IWC  near 
1000  LT  where  SNOE  observes  and  near  2300  LT  where 
SOFIE  observes  is  a  factor  of  3.4.  This  underscores  the 
importance  of  including  tidally  induced  IWC  variations 
when  comparing  data  sets  with  each  other  and  with  model 
results. 


[39]  To  our  knowledge,  Figure  7  shows  results  of  the  first 
PMC  IWC  calculation  over  the  diurnal  cycle  and  has 
important  implications  for  reconciling  two  or  more  sets  of 
PMC  observations  at  discrete  LT.  In  general,  the  IWC  has  a 
maximum  between  0700  and  0800  LT  and  a  minimum 
between  1900  and  2200  LT  with  a  variation  exceeding  a 
factor  of  5  Additional  ground-based  or  satellite  observations 
of  IWC  over  the  diurnal  cycle  would  greatly  help  us  to 
constrain  the  variability  indicated  by  the  simulations. 
Analysis  of  the  Southern  Hemisphere  as  well  as  other  sea¬ 
sons  and  latitudes  would  also  help  to  determine  how  rep¬ 
resentative  these  results  are.  We  now  consider  observations 
from  the  same  June  2007  time  period  at  57°N. 

3.2.  Results  for  Subpolar  Latitudes 

[40]  One  additional  data  set  that  can  be  used  to  validate 
the  NOGAPS-ALPHA/CARMA  simulations  is  from  the 
SHIMMER  instrument  [Stevens  et  al ,  2009].  SHIMMER 
retrievals  do  not  include  IWC  versus  LT  but  SHIMMER 
nonetheless  observed  two  pronounced  occurrence  frequency 
peaks  near  0600  LT  and  1800  LT  during  the  NH  summer 
2007.  We  show  the  calculated  IWC  at  57°N  in  Figure  8a, 
where  we  have  initialized  the  model  with  the  smoke  distri¬ 
bution  from  Bardeen  et  al  [2010]  at  58°N,  included  water 
vapor  freeze  drying  as  in  Figure  7  and  used  the  same  ice 
particle  lifetime  as  in  Figure  6c. 

[41]  Figure  8a  shows  a  strong  IWC  peak  near  0500  LT  in 
good  agreement  with  the  SHIMMER  frequency  peak  near 
0600  LT.  Variations  in  the  assumed  ice  particle  lifetime 
change  the  absolute  IWC  of  the  0500  LT  peak,  but  the 
pronounced  variation  of  the  IWC  calculated  over  the  diurnal 
cycle  does  not  change.  Simulations  that  include  LT  varia¬ 
tions  at  57°N  latitude  and  without  meridional  transport  yield 
no  clouds  at  all.  We  conclude  that  this  early  morning  peak  is 
therefore  due  to  either  southward  meridional  transport  from 
higher  latitudes  where  the  air  is  colder  (Figure  2b)  [e.g., 
Gerding  et  al ,  2007]  or  transient  colder  periods  near  the 
temperature  minimum  not  reflected  in  the  monthly  average. 

[42]  The  second  SHIMMER  peak  near  1800  LT  does  not 
appear  in  the  simulation  of  average  June  2007  conditions 
due  to  the  higher  temperatures  and  subsaturated  air  during 
this  time  of  day.  However,  average  conditions  may  not  be 
representative  for  the  relatively  low  cloud  frequencies 
observed  by  SHIMMER,  which  are  less  than  5%  averaged 
over  all  LT  for  the  northern  2007  summer  [Stevens  et  al , 
2009].  We  herein  consider  the  possibility  that  the  second 
peak  in  Figure  8a  may  be  the  result  of  temperatures  occa¬ 
sionally  falling  below  the  frost  point,  even  though  the 
monthly  average  would  not  reflect  such  cold  temperatures. 
We  represent  this  variability  by  the  standard  deviation  of  the 
NOGAPS-ALPHA  temperatures  at  83  km  altitude  for  the 
month  of  June.  Figure  8b  shows  the  degree  of  supersatura¬ 
tion  (S)  at  83  km  altitude  as  well  as  S  after  subtracting  the 
standard  deviation  from  NOGAPS-ALPHA  temperatures 
during  each  indicated  LT.  The  SHIMMER  PMC  observa¬ 
tions  were  made  between  50°N  and  58°N,  so  we  use  54°N 
as  a  representative  latitude  for  our  analysis  of  supersatura¬ 
tion.  As  one  can  see,  when  the  temperature  falls  below  the 
average  consistent  with  the  standard  deviation  of  the  vari¬ 
ability  the  air  can  often  be  supersaturated  with  values  locally 
peaking  near  three  at  1800  LT.  We  note  that  there  is  equa- 
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Figure  8.  (a)  The  calculated  IWC  at  57°N  for  June  2007 
(black).  Also  shown  for  comparison  is  the  observed  relative 
PMC  frequency  during  NH  2007  from  SHIMMER  [Stevens 
et  al. ,  2009].  The  shaded  region  indicates  where  there  is  not 
enough  sunlight  for  SHIMMER  to  reliably  infer  the  fre¬ 
quency.  (b)  Calculated  degree  of  supersaturation  (S)  over  the 
diurnal  cycle  at  54°N  and  83  km  geometric  altitude  (black). 
Ice  particles  can  exist  where  S  >  1 .  The  blue  curve  indicates 
the  degree  of  supersaturation  after  subtraction  of  the  stan¬ 
dard  deviation  of  the  temperature.  The  SHIMMER  relative 
frequencies  from  Figure  8a  are  again  overplotted  in  red  and 
referenced  to  the  right  hand  axis. 


torward  flow  until  about  1500  LT  (Figure  2b)  with  a  merid¬ 
ional  wind  reversal  near  this  time  so  that  a  peak  in  cloud 
frequency  might  also  be  expected  near  this  time  from  equa- 
torward  transport.  The  NOGAPS-ALPHA  analysis  indicates 
therefore  that  conditions  can  often  be  favorable  for  the 
appearance  of  clouds  at  1800  LT  even  though  the  monthly 
average  result  in  Figure  8a  does  not  show  this. 

4.  Discussion 

[43]  Our  study  on  tidally  induced  variations  of  PMC  IWC 
for  June  2007  approaches  the  problem  in  three  new  ways. 
First,  we  have  chosen  to  simulate  the  IWC  rather  than  cloud 
brightness  or  frequency  for  four  reasons:  (1)  IWC  is  a  direct 


measure  of  how  much  ice  is  present,  (2)  it  is  a  vertically 
integrated  quantity  and  therefore  does  not  require  simulation 
of  a  narrow  ice  layer  peak,  (3)  it  is  primarily  sensitive  to 
only  temperature  and  water  vapor  [Rapp  and  Thomas , 
2006],  and  (4)  it  is  measured  directly  by  SOFIE  on  the 
AIM  satellite.  Second,  we  have  used  meteorological  analy¬ 
sis  fields  from  a  global  numerical  weather  prediction  system 
that  has  recently  been  extended  to  ~90  km  altitude.  This 
allows  us  to  specify  temperature,  water  vapor,  and  winds 
over  the  diurnal  cycle  for  a  specific  time  period  during  which 
we  have  PMC  data.  Third,  we  drive  the  one-dimensional 
version  of  CARMA  with  time  varying  NOGAPS-ALPHA 
input  to  simulate  PMC  altitude  variations  and  IWC  over  the 
diurnal  cycle.  Specification  of  assimilated  ambient  condi¬ 
tions  with  latitude  and  LT  more  realistically  represents  the 
ice  particle’s  history  than  either  a  free  running  climate 
model  or  a  microphysical  model  with  a  constant  temperature 
profile  as  is  sometimes  done.  We  assume  that  only  cosmic 
smoke  particles  serve  as  nucleation  sites  for  the  clouds  and 
that  clouds  formed  in  any  other  way  contribute  negligibly  to 
the  total  IWC  at  any  given  LT. 

[44]  One  limitation  of  this  study  is  that  it  uses  assimilated 
data  from  only  the  month  of  June  2007,  which  may  not 
represent  typical  June  conditions.  To  address  this  limitation, 
we  have  continued  the  NOGAPS-ALPHA  assimilation 
configuration  described  in  the  work  of  Eckermann  et  al. 
[2009]  through  the  2009  PMC  season  and  assembled  the 
temperatures  for  three  consecutive  Junes  on  a  zonally 
averaged  LT  grid.  Figure  9  shows  the  variation  of  temper¬ 
ature  over  the  diurnal  cycle  at  83  km  geometric  altitude  for 
June  2007,  2008,  and  2009  at  three  different  latitudes:  81°N 
(Figure  9a),  69°N  (Figure  9b),  and  57°N  (Figure  9c).  The 
consistency  of  the  assimilated  temperatures  from  one  year  to 
the  next  is  remarkable  and  on  average  the  temperature  is 
about  10  K  lower  at  81°N  compared  to  57°N. 

[45]  With  respect  to  the  temperatures  over  the  diurnal 
cycle,  we  draw  several  important  conclusions  from  Figure  9: 
(1)  the  migrating  diurnal  tide  dominates  LT  temperature 
variability  throughout  the  Arctic  mesosphere  in  June  from 
2007  to  2009;  (2)  the  average  amplitude  of  the  tidal  tem¬ 
perature  oscillation  is  repeatable  at  any  given  latitude  for 
this  time  period  and  the  time  of  the  maximum  varies 
between  1400  and  2000  LT  at  all  latitudes;  (3)  the  amplitude 
becomes  progressively  larger  at  lower  latitudes  with  an 
amplitude  of  ~2  K  at  81°N  to  ~5  K  at  57°N;  and  (4)  the 
temperature  variability  at  a  given  LT  also  becomes  pro¬ 
gressively  larger  at  lower  latitudes  with  a  standard  deviation 
of  3-6  K  at  81°N  and  5-8  K  at  57°N.  We  conclude  from 
Figure  9  that  the  diurnal  variation  of  temperatures  for  our 
June  2007  simulation  is  representative  for  solar  minimum 
conditions  at  all  latitudes.  Since  we  find  that  the  diurnal 
variation  of  temperatures  primarily  controls  the  IWC  vari¬ 
ation,  we  also  conclude  that  the  calculated  IWC  over  the 
diurnal  cycle  shown  in  Figure  7  is  also  representative  of 
solar  minimum  conditions. 

[46]  Ultimately,  our  analysis  suggests  that  our  work  may 
have  the  greatest  utility  in  interpreting  data  from  different 
satellites  which  are  in  Sun- synchronous  orbits  fixed  in  LT. 
As  we  noted  in  the  introduction,  the  SBUV  series  of 
instruments  has  viewed  PMCs  for  over  30  years,  which  is 
the  longest  PMC  time  series  available.  These  nadir-viewing 
observations  are  typically  reported  as  albedo  (unitless)  but 
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Figure  9.  (a)  The  variation  of  temperature  over  the  diurnal 
cycle  for  June  2007-2009  at  83  km  geometric  altitude  and 
81°N.  The  vertical  bars  represent  the  standard  deviation  of 
the  assimilated  temperatures  within  the  monthly  averages, 
(b)  Same  as  in  Figure  9a,  except  for  69°N.  (c)  Same  as  in 
Figure  9a,  except  for  57°N. 


are  actually  directional  albedo  (sr_1)  at  252  nm,  with  mul- 
tidecadal  trends  of  less  than  1%/yr.  Figure  10  shows  the 
variation  of  directional  albedo  at  252  nm  for  each  hour  of 
LT  using  the  same  CARMA  size  distributions  as  the  IWC 
calculations  in  Figure  6c.  We  calculate  the  albedo  at  a 
variety  of  relevant  SBUV  solar  scattering  angles  (SSAs) 
from  90°  to  150°.  Away  from  the  terminator  ( 0  >  90°),  the 
average  calculated  albedo  in  Figure  10  is  between  0.2  and 
3  x  10~5  sr-1,  consistent  with  the  albedo  typically  reported 
for  SBUV  of  about  1  x  10~5  sr-1  [e.g.,  DeLand  et  al .,  2006]. 

[47]  Figure  10  illustrates  that  the  calculated  albedo  vari¬ 
ation  over  the  diurnal  cycle  is  similar  to  the  IWC  variation, 
with  a  peak  between  0600  and  0800  LT  and  a  minimum 
between  1900  and  2100  LT.  The  largest  variations  over  the 
diurnal  cycle  reported  heretofore  are  in  cloud  occurrence 
frequency,  which  are  implicitly  included  in  our  averaged 
model  results.  Large  variations  in  cloud  frequency  over  the 
diurnal  cycle  with  a  prominent  peak  in  the  morning  were 


reported  in  both  ground-based  and  SBUV  data  sets  [Fiedler 
et  al .,  2005;  Shettle  et  al. ,  2009]. 

[48]  It  is  well  known  that  the  mid-UV  PMC  albedo  is  a 
function  of  the  SSA  [van  de  Hulst ,  1981],  and  we  note  the 
SBUV  SSA  is  the  supplement  of  the  solar  zenith  angle.  The 
curves  show  that  the  albedo  is  significantly  larger  at  the 
terminator  (6  =  90°)  than  for  lower  solar  zenith  angles  at  a 
given  LT.  This  means  that  a  comparison  of  SBUV  albedo 
and/or  frequency  observations  not  only  requires  consider¬ 
ation  of  the  LT  but  also  the  SSA,  particularly  when  the 
observations  are  near  local  sunrise  and  sunset.  Importantly, 
the  calculated  albedo  (or  IWC)  maximizes  in  the  morning 
which  corresponds  to  the  average  LT  of  the  SBUV  observa¬ 
tions  between  1979  and  1990  after  which  the  observations 
are  more  uniformly  spread  over  the  diurnal  cycle  [DeLand 
et  al .,  2007].  If  the  June  2007  results  in  Figure  10  are 
representative  of  all  northern  PMC  seasons  within  the 
multidecadal  SBUV  data  set,  this  implies  that  inferred  long¬ 
term  PMCs  trends  of  albedo  or  IWC  would  increase  fol¬ 
lowing  adjustment  due  to  diurnal  variations. 


5.  Summary 

[49]  In  summary,  we  have  produced  the  first  quantitative 
simulation  of  IWC  over  the  diurnal  cycle  at  57°N  and  69°N 
using  assimilated  satellite  data  from  June  2007.  The  simu¬ 
lation  at  69°N  is  constrained  by  concurrent  observations 
from  SOFIE,  which  measures  IWC  directly  near  23  LT.  The 
IWC  has  a  diurnal  variation  that  is  controlled  by  the 
migrating  diurnal  temperature  tide  and  the  peak  IWC  is 
between  0500  and  0800  LT  at  both  57°N  and  69°N. 

[50]  The  relative  variation  in  IWC  over  the  diurnal  cycle  is 
larger  at  the  equatorward  fringe  of  the  PMC  region,  signif¬ 
icantly  complicating  the  interpretation  of  multidecadal 
trends  reported  at  these  latitudes  [DeLand  et  al .,  2007]. 
Even  at  69°N,  however,  the  IWC  or  mid-UV  albedo  varies 


Figure  10.  Calculated  PMC  directional  albedos  for  obser¬ 
vations  at  252  nm  at  a  variety  of  SSAs  ( 0 )  relevant  to  the 
SBUV  data  set.  The  calculations  assume  spherical  particles 
and  use  the  same  particle  size  distributions  as  the  IWC 
calculations  in  Figure  6c.  For  the  nadir- viewing  geome¬ 
try  of  SBUV,  the  solar  zenith  angle  is  equal  to  1 8O°-0  so 
that  0  =  90°  is  representative  of  observations  at  the  Earth’s 
terminator. 
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by  at  least  a  factor  of  5,  with  a  peak  between  0700  and 
0800  LT  and  a  minimum  between  1900  and  2200  LT. 
Careful  consideration  of  not  only  the  LT  but  also  the  latitude 
for  any  given  set  of  observations  must  therefore  precede  any 
IWC  trend  analysis. 

[51]  Our  description  of  IWC  over  the  diurnal  cycle 
currently  lacks  experimental  validation.  Ground-based 
observations  of  the  PMC  frequency  at  69°N  show  a  factor  of 
4  variation  with  a  peak  near  1  LT  [Fiedler  et  al ,  2005]. 
Limited  satellite  observations  of  PMC  frequency  are  con¬ 
sistent  with  a  factor  of  3  variation  during  the  morning  with  a 
maximum  near  2  LT  [Shettle  et  al ,  2009].  Modeling  work 
by  Jensen  et  al  [1989]  showed  that  the  PMC  brightness 
varied  by  up  to  a  factor  of  7  over  the  diurnal  cycle 
depending  on  the  amplitude  of  the  temperature  variations, 
with  a  maximum  near  23  LT  that  follows  a  temperature 
minimum  near  20  LT.  Although  none  of  these  studies  report 
the  IWC  variation  over  the  diurnal  cycle,  taken  together  they 
illustrate  that  variations  of  related  quantities  can  be  between 
a  factor  of  3  and  7.  Our  derived  variation  of  at  least  a  factor 
of  5  in  IWC  at  69 °N  is  therefore  not  unreasonable  in  the 
context  of  these  previous  studies.  Our  inferred  maximum 
between  0700  and  0800  LT  is,  however,  between  5  and  9  h 
later  than  these  results  would  indicate.  Direct  temperature 
observations  near  83  km  altitude  over  the  diurnal  cycle  in 
addition  to  any  IWC  observations  would  provide  a  useful 
comparison  to  the  results  presented  here. 
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